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SCME CONICAL AND QUASI -CONICAL FLOATS IN LINEARIZED 
SUPERSONIC -WING THEORY 
By Herbert S. Ribner 


SUMMARY 

A number of conical and quasi -conical linearized supersonic 
flows have been derived „ These flows may be applied in lift- 
cancellation techniques in the determination of wing-lift dis- 
turbances that arise at subsonic trailing edges. Some of the 
results are applied to damping in roll and pitch in another paper. 

Two methods of an al ysis have been employed, both involving 
integral equations. One is a development of the membrane method 
of Eward; the other is an analog of the superposition method of 
Schlichting. In both cases, source distributions are used. The 
mathematical equivalence of the first method and a new doublet - 
distribution method of Goodman and Mirels is shown. 


INTRODUCTION 

The analysis of tie flow over a sweptback wing (reference 1) 
makes use of Lagerstrom's concept of lift cancellation (reference 2). 
The starting- point is a delta wing of infinite chord. Suitable flows 
are superimposed that cancel the lift outboard of and behind a 
certain bo unda ry. This boundary is chosen to constitute the tips 
and the trailing edge of the sweptback wing. 

The cancellation flows for the wing tip modify the lift in the 
tip region; they are not of concern herein. The cancellation flows 
for the trailing edge modify the lift in a region ahead of the trail- 
ing edge if the edge is "subsonic"(that is, if the component stream 
velocity normal to the edge is subsonic). Now the lift distribution 
of a delta wing at an angle of attack is substantially flat in the 
center. (See fig. 1.) Thus, an approximate cancellation behind the 
trailing edge would be afforded by a constant lift such as flow I in 
the figure. (This idea was originated in reference 1.) Correspond- 
ing approx imat e cancellation flows (ill and IV) for rolling and 
pitching motions are also shown in the figure. The derivation of 
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these three flows is the main object of this report; the work was 
initiated at the NACA Langley laboratory and completed at the NACA 
Lewis laboratory. 

The analysis leads to integral equations, some with several 
apparent solutions. The choice of the physically correct solution 
in each case is determined by comparison with the solution for a 
closely related flow that is more easily calculated. Thus, 
altogether, derivations are given for eight flow distributions. 
Four of these flows are conical: the velocity components and pres- 
sure are' constant along any ray from the vertex. The remaining 
four are quasi -conical: the velocity components and pressure are 
proportional to the distance along any ray. 

The practical application of several of these cancellation 
flows in the evaluation of damping in roll and pitch is made in 
reference 3 . 


CANCELLATION FLOWS 

Each of the eight flows discussed can be considered to repre- 
sent a slightly cambered lifting surface lying essentially in the 
2 = 0 plane. (All symbols are defined in appendix A.) Thus, the 
u and v velocities are antisymmetric with respect to the z = 0 plane, 
and the w velocity is symmetric. With these reservations, the 
boundary conditions can be specified for the upper surface only. 

(This simplification is employed throughout the report.) 

The flows contemplated for the partial cancellation of the 
lift behind the trailing edge are schematically designated in 
figure 2 for the several cases. A particular cancellation flow 
(any one of I to IV, fig. 2) is superposed on the sweptback wing 
so that the w = 0 regions lie ahead of the trailing edge and the 
shaded region lies behind. The lift, and hence the u velocity, is 
specified in the shaded region, and this specification determines 
the lift cancellation. (See fig. 1.) The specification w = 0 for 
the part -of the flow overlapping the wing ensures that the 
resultant flow shall be tangent to the surface of the wing; the 
basic delta-wing flow already provides the correct value of upwash w 
on the wing (for example, -dV for angle of attack) and the cancel- 
lation flow must therefore add none . 

The u velocity of the cancellation flow in the regions w = 0 
is algebraically additive to that of the basic delta flow there. 

ThiB additive u velocity corresponds to a pressure disturbance caused 


1317 



NACA IN 2147 


3 


"by the wing trailing edge. The central problem of this report is 
the determination of the u velocity for each of the cancellation 
flows I to IT from the boundary conditions stated in the preced- 
ing paragraph. 

In the integral -equation method of solution as applied herein, 
several solutions may be found for certain of the integral equ- 
ations. Unique solutions can be obtained in each case,’ however, 
for a flow that is closely related to the desired flow. The 
correct solution for the desired flow may be identified by a com- 
parison of the nature and the location of singularities in the two 
flows. The identification is made more convincing by a quanti- 
tative comparison of the u velocity over the region of interest in 
addition to a consideration of the singularities . 

Thus in figure 2 flows I to IV are the desired flcwB and flows 
I' to IT* are the respective related flows. The boundary conditions 
for the desired and related flows differ only in the presence or 
absenoe of a left-hand w = 0 region. Because of the relative remote- 
ness of the left-hand region, conditions in thiB region may be 
expected to affect but slightly the u velocity induced in the right- 
hand w = 0 region, so long as that region is relatively narrow. With 
this limitation the u velocity in the right-hand w = 0 region cal- 
culated for flow I', then, is presumed to be a good approximation to 
that for flow I, and similarly for the other pairs of flows. 

It will be convenient to obtain first the solutions for the 
related flows I* to TV' so that they will be available to aid in 
Identifying the correct solutions for flows I to IT, respectively. 
Flows I* to IT' require only a simple inversion of an Abel -type 
Integral equation. Flows I to IT require a more elaborate procedure. 


Flow I’ 

(For Angle of Attack) 

In flow I' (fig. 3) the w velocity (or surface slope) is 
prescribed in region A, but not the u velocity (or lift). The 
u velooity is prescribed in region B, but not the w velocity. If 
the w velocity were known in region B, this wing would be of pre- 
scribed c amb er - According to Puckett (reference 4), the wing 
could be represented by a distribution of sources in proportion to 
the local value of w. The unknown u velocity in region A could 
then be obtained by simple integration and differentiation. The 
source representation can still be used, however, even though the 
distribution of w in region B is unknown. In this case the pre- 
scribed conditon on u in region B will give rise to a soluble 
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integral equation for the unknown v distribution.' This procedure 
is a development of Eward's method (reference 5). The details 
follow. 


The surface u velocity is given in reference 5, equation (3a). 
In slightly modified form, and in the present notation, it is 


( 1 ) 



The oblique coordinates of reference 5 (Mach coordinates) are employed. 
(See fig. 3.) Unsubscripted values refer to the field point, sub- 
scripted values to the source point. The transformations relating 
(x,y) and (£,q) are 


X = m £ 

Py = | (ti-£) Ti 

The surface integral is taken over the area S bounded by 
the forward Mach cone from the point (!,tj). The line integral is 
taken along any lines within S or bounding S across which the 
w velocity experiences a Jump Aw. The | -axis and the line q = T | 
(fig. 3) might be such lines. Consider, however, the requirement 
of general similarity with flow I. The imposition of the Kutta 
condition at the trailing edge of the wing in figure 1 requires the 
continuity of w across the boundaries q = T | and | = TT) of the 
shaded region of flow I, as sketched in figure 2 . For flow I * , w 
must s imilar ly be continuous along the right-hand edge of the shaded 
region (q =t|). Also, the continuity of w across the i -axis 
(Mach line) can be shown to be a consequence of the finite lift along 
that line. The line integral in equation (l) accordingly vanishes 
for flow I ' . 

Let the point (t,q) be looated at P (fig. 3); P repre- 
sents an arbitary point in region B. Then the forward Mach cone 
inoludes the shaded area in region B and an adjacent unshaded area 
in region A. The quantity Sw/cis^ =0 in region A; therefore, 

only the shaded area contributes to the surface integral. 

Equation (1) may thus be written 


- (x - Pr) 
” C»fir) 


*\ 


J 


( 2 ) 
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The unknown inner integral, multiplied hy l/itM, may he designated 
which gives 


T] 

J(^ T h)^- T l-i 

u(£,q) = - — ~ — region B (4) 

nJ o 


Equation (4) is an Ahel-type integral equation for the unknown 
function J(|,Tj-j_)» solution is given (reference 6) hy 


J(^l) - "7 


<5% 


>7 ll 


0 


u(£>n)foi 


(5) 


(The more general treatment of Ahei’s integral equation in refer- 
ence 7 implies that equation (5) will provide singular as well as 
continuous solutions. The only Important restriction is that 
u(£,q) must he such that the integral on the right-hand side of 
equation. (5) is continuous. In reference 6, however, certain 
additional restrictions axe placed on u(!,q) to limit equation (5) 
to continuous solutions. In the present report, these additional 
restrictions are disregarded and singular solutions may he expected, 
as in equation (6).) 


In the present case the function u(I,t]) 
(fig. 3). The solution for J is therefore 


is a constant Uq 


J(Ml) 



( 6 ) 


Wot let the general point (£,t)) he located in region A 
at Q (fig. 3). The only change in the expression for u 
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(equation. (3)) is in the upper limit for t^. This upper limit is 

now at the intersection of the line PQ with the line ij = t| . 
Accordingly, 

't£ 

u(£,q) = - | — region A (7) 


0 




The -value of J already obtained in equation (6) applies here as 
well as in equation (4) . Upon making the substitution, 

dT b. 

nni -^ 2 


= 5^2 sin -1 \[lL region A (8) 

it V T1 


u 0 

u(S,n) = T 






It will he convenient to' reexpress the result of equation (8) 
in terms of the conical coordinate O - 3y/x. Along the line 
q = t| , the coordinate <J is given the value n. The transformation 


I 






^ =• Hr x < 1+a > ) 



(9) 


( 10 ) 


Equation (10) gives the hitherto unknown u velocity in region A 
of figure 3. The w velocity in region B is still unknown. This 
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v velocity -will not "be needed in the present investigation, "but 
the method of solution is given in appendix B as a matter of 
interest . 


Flow H’ 

(Special) 

The development leading to equation (5) in the treatment of 
flow I* leaves the "boundary condition on u as yet unspecified. 
The development and' the equation are thus sufficiently general to 
he applied to any of the flows I' to 17'. 

The present case is complicated hy discontinuity in the 
u velocity across the x-axis (fig. 4). The point P = (£,t)) may 
lie in either of the regions B^ or B£ . The corresponding 
solutions for J ( from equation (5) will he different; c al l 
them Jq(l > and respectively. The conditions on 

u(|,n) are 

u(|,tj) = -u 0 0 < ti < | region B x 

u(t,T]) = vlq i < n < T| region Bg 


Then 


^ 








pni 


-UQdTl 

Yt 1-n 


u 0 dT l 
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The results are 






% u 0/ 1 2 
J 2 (£ ’" l) = « ( ^ ' Vv 7 ! 


Equation (7) is, like equation (5), applicable to any of the 
flows I' to IF’ . For flow II ' the appropriate function Jd,^) 
is JiCljTh) in the range 0 < n < t and in the 

range £< tj < t£. Therefore 


= - 


J 1 (t,Ti 1 )dq 1 


J 2 ( ^l )dT h. 


+ 2 


In^i-v 




2u 0 _l (2T-l)£-q _j_-l . T l region A 

= rt - cos -*■ — -=tE — *■ - Bin \ ~ 

It T]“ £ V ^ 


In terms of the conical coordinate O = Py/x - this equation is, 
after the transformation equations (9) are applied. 


1 2 n-rtl «l . sin -l,HS 


— “ " 003 0 (i-n) ' \|(i-n)(l+OJ 
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Equation (12) gives the u velocity in the v = 0 sector 

( n < a < !) of flow II* (fig* 1), which is zone A of figure 4. 


plow hi' 

The appropriate value of u to he substituted in equation (5) 
is, from figure 2, 


u = icy 


or 


u = I h-V 


so that 


J(l,T)l) = - jtM 


^1 


.(HtDA 1 ! 

V 11 ! 71 


JL 

* M ^ T 


(13) 


Insertion of this value into equation (7) gives 

prl 




(S^-tJdT^ 

^lC^i-ni) 


(T} > t|) 


(14) 


which integrates to 


u(|,Tl) = 


2_jc 

JtM 


(t)-£) sin 


-1 


N 



t£(tj-t! ) 
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In terms of the conical coordinate G = Py/x this equation is, 
after the transformation equations (9) are applied. 


u(g,x) 



V 2 (1+n) (1 -g) (q-n) 
1-n 


(15) 


Equation (15) gives the u velocity in the w = 0 sector 
(n < O < l) of flow UI ' of figure 2. 


Flow 17* 

Refer to figure 2 for the value of u to he substituted 
in equation (5) : 


u = Kx 


“ f- (*•» 


Thus 




j(t,Tl 1 ) 


Kp d 

^5^ 


(m-PdTi 

i\fn 


Kp 2r ii + ^ 
* M 


(16) 


Substitution of this value in equation (7) gives 

r> t| 


u(|,ti) = 


Kp 

jtM 


(„ > r£) 


The integration yields 

X 2*£ 

U(|,ll) = rtM 


(*n+l) si n ' 


' 1 \|t • ^ 


1-Tt) 
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In terms of the conical coordinate a = Py/x this equation is, 
after the transformation equations (9) are applied. 


u(o,x) 


2 8in“^ / 
2 sin 

+n) (l-a) M 2 (l+n) (l-a) (cr-n) 


-n) (1+a) 1-n 


(17) 


This equation gives the.u velocity in the w = 0 sector 
(n < O < 1) of flow IT' of figure 1. 


Flow I 

(For Angle of Attack) 

Flows I to IT differ frcm flows I* to IT', respectively, in 
the specification of a left-hand sector w = 0 to match the right- 
hand sector (fig. 2). The solution can no longer he obtained by 
the simple inversion of an Abel -type integral equation. Resort is 
therefore made to a different integral -equation formulation in 
which advantage is taken of the conical nature of the flow. The 
method is an analog of the superposition method of Schlichting 
(reference 8) . 


According to the considerations developed in the discussion of 
flow I', flow I may be represented by a suitable source distribution. 
At every point the required source strength is proportional to the 
local value of w. Thus the source strength is zero in the two 
outside sectors. (See fig. 2.) In the central sector w is 
u nkn own, but there is a condition on the u velocity. The source 
distribution there must be so chosen that this condition on the 
induced u velocity is met. The formulation of this condition gives 
rise to an integral equation for the unknown source distribution. 

A convenient form of integral equation is obtained as follows i The 
source distribution over the region u = uq is considered to be 

obtained by superposing uniform sectors of infinitesimal strength 
and different lateral extent Section A -A of figure 5 illus- 

trates such a distribution. The strength of an elementary sector is 
given by 


by 

5ci-[ 



Aw = 


(a x > o) 
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Each elementary source -sheet sector induces a certain infinitesi- 
mal u velocity . and the total u velocity is given "by the integral. 
This integral must he equated to the prescribed u velocity shown 
in figure 5. The result is the desired integral equation for the 
unknown distribution strength Sw/Bc^ . 


The infinitesimal u velocity induced by an elementary source - 
sheet sector extending from -O^ to can be obtained from 

reference 4 (equation (31) ) . In the present notation it is 


Au(o) = 



i-c^a 

a "°i| 


+ cosh 


-1 



(18) 


where an absolute value sign has been added to the denominators of 

the cosh -1 terms. With the absolute value sign, equation (18) 
applies both for |a| < and for |a| > whereas in reference 4 

the two cases are separately covered in equations (31) and (33) } 
respectively. (The same expression results from the addition of two 
oppositely swept line source's of the acceleration potential. See 
reference 9, equation (12).) 

The total u velocity induced by the superposed source -sheet 
sectors is 


n 




Sw 

^1 


Ol 




cosh 




-i 

la -a, 


+ cosh" 


1+OiL 

"F?l[ 


\ 

7 


(19) 


If u(a) is put equal to the specified value Uq in the range 
|a| < n, then equation (19) represents an integral equation for 
the derivative dwta^/Sa^ of the unknown source strength. This 


equation as it stands is too complex to be useful. Great simpli- 
fication results, however, upon differentiating both sides with 
respect to o: 


nn 


du 

3a 


2a 


up 




isr °i do i 




2 -o 2 


( 20 ) 
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(The Cauchy principal value is to he taken for the integral.) In 
the range (a) < n, du/Sa = 0. The following less formidable 
integral equation is thus obtained: 


'n 


dy 

3oJ CT i 


da x = 0 |CT| < n 


( 21 ) 


2 2 

The substitution tp = Qp , t = O reduces equation (21) to a 
special case of 


Pt 


^a 


t-^-t 


= g(t) 


(21a) 


The integral equation (21a) is well known from incompressible thin- 
airfoil -section theory, but the usual inversions yield only the 
trivial solution f(tj) = 0 for g(t) =0. A nontrivial solution 

of equation (21) for dw/Ba^ has been suggested by C_. E. Brown of 
the BACA Langley laboratory. Put 


C? = ~ (1 - cos 8) 

= y- (! “ cos e ) 




> ( 22 ) 


ST " cf(e) 


(C = unknown constant) 


J 


Then equation (21) becomes 


f(9) sin 9 d9 _ q 
cos 9 - cos 8 


( 23 ) 


0 
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Now according to reference 10 

r* 


&e 

cos 9 - cos 5 


= 0 


whence a solution of equation (23) is 


or 


f(0) = osc 9 




dw _ 

'S5[~ 


= C osc 9 


J 


(24) 


The question of uniqueness is deferred until later; in the meantime 
equation (24) is considered as the physically correct solution. 

Equation (24) carries the solution for far enough to 

enable the determination of u. Thus, equation (20) is still valid 
in the range n < |crl <_1. In this range Su/ckJ ^ C, and 
equation (20) may serve for the evaluation of u. The transformation 
from a 2 to 9 (equations (22)) is again convenient. (The trans- 
formation from a to 6 is not used as it leads to cos 5 > 1.) 

The result is 




du _ CO 




0 


f(8) sin 9 dO 
2 

— 1 + cos 9 

n 2 


n <|a| < 1 


and with f(0) = esc 9 


nit 


du _ 

$5" 


Ca 


itp^l-c? 


2a 2 . 


d£ 

1 + cos 9 


n<|a|< 1 (25) 


JO 
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This integral is evaluated in reference 11 (table 64, equation (12) ) . 
The result is 


3u 

35 = 


Cn^O 

26 ct|^(1-c^) (a 2 -n 2 ) 


n < |CT| < 1 


(26) 


For purposes of integration, c may be restricted to the positive 
range so that cr/|a| =1; the resultant integral -will be applied 
for both the positive and the negative value of O because of the 
symmetry of the flow in O . The elliptic -function substitutions of 
appendix C are helpful in this and other more difficult integrations 
of the same nature in this report. The integral of equation (26) 
between the limits a and 1 is found to be 


u(a) = 


F(<P,k) 


where F is the incomplete elliptic integral of the first kind 
with modulus 


and amplitude 



cp 



The constant C is established by the condition that u(n) = Uq 
bo that finally 


u(a) = Uq 


F(cp,k) 



(27) 


Equation (27) is a solution for the u velocity in the out- 
board sectors n <_ O < 1 of flow I in figures 2 and 5. Equation (27) 
is not a unique solution of the integral equation (21) , and it remains 
to be shown that it is the physically correct solution for flow I. 

The singularities in 3u/3a, equation (26), may be compared with 
those in the corresponding result for flow X', obtained by differenti- 
ating equation (10) . In both cases half -order singularities are 
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found at o - n+ -and o = 1 (that is, |du/3a| approaches infinity 

as (a-n)~^/ 2 when O approaches n from above (designated hy 
rn-) and there is a similar behavior as o approaches 1) . The 
quantitative comparison of equation (27) with the corresponding 
solution, equation (10), for the right-hand sector of flow I' 

(figs. 2 and 3) is even more convincing. The two solutions are 
plotted together in figure 6 for n = 0.707. The agreement is so 
close that the two sets of points seem to define a single curve, the 
uppermost curve of figure 6. (Other calculations for the point 
defined hy <p = 60°, plotted in fig. 7, confirm the expectation 
that this close agreement becomes progressively impaired as n is 
chosen smaller and smaller. The error is within 5 percent down to 
n = 0.32 .) 

Equation (27) is the correct solution for the "symmetrical 
wake correction" of. reference 1. This result has been incorporated 
in an erratum sheet thereto. 


Flow U 
(Special) 


Flow II has more academic than practical significance . It 
could be used in the determination of the loading cn .special ailerons 
to cancel (approximately) the loading in the wake of a sweptback wing 
with the ailerons deflected. These ailerons would be full span, or 
located inboard if part span, with a vertical fence at their juncture 
to isolate the two ahead of the region of trailing -edge disturbance. 


Comparison of the- specifications of flow H and flow I in 
figure 2 shows that II is antisymmetric whereas I iB symmetric. It 
is therefore necessary to use antisymmetric source -sheet sectors in 
the superposition process, rather than the symmetric sectors used 
for flow I. The change is effected by simply changing the sign of 

the second cosh -1 term in equations (18) and (19) . The equation 
corresponding to equation (20) then becomes 


r>n 


3a 


rtP^l-a 2 


a, 2 da 1 
da x 1 - 1 

o^-o 2 


( 28 ) 


Jo 
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(The Cauchy principal part of the integral is to he taken.) The 
Integral equation that results on setting Su/Sa = 0 is thus 


r a 


Sw 

*5 


On C dOl 


a-^-o 2 ' 


= o |a[ < n 


(29) 


in place of equation (21). Equation (29) is of the same form as 
equation (21) : dw/3o-]_ is replaced by as the unkn own 

function of 0^* According to equation (24), a solution of 
equation (29) therefore is 


chr 

35[ 


(Jn = C esc 9 


(30) 


in terms of the function 0 defined in° equations (22). 

For the determination of u, the solution equation (30) is 

2 

substituted back into equation (28) with CJ^ 2 = — (1 - cos 0) 

2 

as in equation (22). There is obtained 




Su 

cKJ ~ " 



o' 


0 


d0 


i£ - i + 


008 0 


n < \o\ < 1 


This expression for chi/Ba differs from equation (25) only in 
the lack of a factor O outside the integral sign. Comparison 
with equation (26) therefore shows that the result of the inte- 
gration must be . 


Gn 


2P|a| ^(i-o2)(a2- n 2) 


Su 

S5 = 


n < |o| < 1 


(31) 



o 
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This function is readily integrated with the aid of the elliptic- 
function substitutions of appendix C. (The absolute value sign is 
temporarily ignored.) The integral from a to 1 can be expressed 
in the form 

u(ct) . cob- 1 U+aV.-jin 2 . n < |o| < 1 

« M (1-rfV - - 


where the sign factor a/ |a| has been appended in order to provide 
the antisymmetry of u in ±0 indicated by equation (31) . The 
constant C is established by the condition u(n) = Uq, so that 
finally ■ 

u(CT) = ^2 -£_ cos" 1 ,(l+n 2 )Q 2 -2n 2 n < |a| < 1 (32) 

31 | a I (l-n 2 )a2 


Equation (32) is the solution for the u velocity in the w = 0 
sectors of flow H (fig. 2). As before, the solution is confirmed 
by comparison with the solution (equation (12)) for the related 
flow H' of figure 2, with attention to the singularities in du/dcf. 
The two solutions are plotted together in figure 6 for n =» 0.707. 
Again, the agreement is such that the two setB of points seem to 
define a single curve, the second from the top in the figure . 


Flow HI 

The boundary conditions (fig. 2) show that flow 111 cannot 
be conical like flows I and XI, but must be quasi-conical. That 
is, the velocity components are not of the form u = F(y/x) and 
w = G(y/x), but are rather of the form u = xF(y/x) and 
w = xG(y/x) . The flow is, however, antisymmetric like flow H. 

A flow of this type can be built up (compare treatment of flows I 
anfl tt ) by superposing elementary source -sheet sectors that have a 
strength variation proportional to y or, what amounts to the 
same thing, to xcr. 

An example of how an arbitrary antisymmetric quasi-conical 
w velocity distribution of the form 


w(o L ) = xy-jffcF-].) 


(33) 
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can be built up by such sheets Is shown in figure 8. (Compare 
fig. 5, section A-A.) The strength of an elementary sheet is 
characterized by the strength parameter 




The incremental u velocity induced by an elementary sheet for 
which 


w = (constant) Py^_ 

= (Af) ro 1 

can be derived by standard methods (for example, reference 5) . 
The expression is 


Au 


df 

<302 


3 


• trCa^) 


(34) 


where 


U(0,0 X ) 


a l 2 


*(1-0^) 


2<W 


- *i l+O-i O -» 

(1+0,0) cosh =—r - (1-0, 0 ) COBh" 

|o 1+ CT| 


1-OjO 
02 -0 


The total value of u contributed by the superposition of all the 
contemplated sheets is given by 

° on 

u(x,0) = - "p - 20~ 31(0,02) ^1 (35) 

Jo 


Now u(x,o) is specified over the shaded area (|o| < n) of 
flow IH in figure 2. Thus, equation (35) represents an integral 
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equation for the unknown function df(o-j_)/cl Gq. The kernel 
TJ(G,G^) of this integral equation is quite complex. The second 
derivative of u(o,Oq) with respect to G has, however, the 
much simpler form 


d 2 U(a,0 1 ) _ 400} 3 

3a 2 “ \(a 1 2 -o 2 ) 2 ’i/W 

* 


This simple function would then he the kernel of the integral 
equation obtained by differentiating equation (35), provided that 
the right-hand side of equation (35) could be differentiated under 
the integral sign without regard for the singularity at Gq = g. 

It is clear that such a procedure' would be invalid in this case 
because for |G| < n the kernel would have a second-order pole 
within the range of integration. 


H. Mirels of the Lewis laboratory has, however, pointed out 
that if the differentiation were conducted properly the result 
might be written in the form 


« - **** (36 

So 2 P J d0 l *(0,2- o 2 ) 2 ‘Jl-if 

r 

The symbol signifies that the indefinite integral evaluated 

Jo 

at Gq= 0 is to be subtracted from the indefinite integral 

evaluated at G 1 = n. Defined in this way, the integration can 

be carried out without’ -difficulty, inasmuch as the integrand is 
regular in the neighborhood of the limits . 


The formulation of equation (36) is a special case of a 
theorem that may be briefly stated as follows: If the function 

f(x,£) has an "integrable singularity" within the interval of 
integration, then, under rather general circumstances, the equation 
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■'ll 

-^f(ac,S)a| (D4) 

Sx n 

a 


1b correct. A general discussion and proof of this theorem is 
provided in appendix D by F. E. Moore of the Levis laboratory. 

The prescribed condition on u in the range jo| < n is 
(flow HI, fig. 1): 


a n 

ax n 


f(x,|)d| 


or 


u = Ky ^ 


Kxg 
= P 

> 


^u 


0 


J 


|a| < n 


Application of this condition to equations (35) and (36) yields 
the two alternate Integral equations for df(a 1 )/da 1 




KO = 


df 

dcr. 


U(a,a 1 )da 


} M < n 


r n 


0 = J 

Jo 


df 

dt l (Oj 2 -#) 2 


J 


(37) 


( 38 ) 


Attention will be centered on solving the simpler of the two 
equations, (38). 



22 


mCA TN 2147 


It is convenient to .-write equation (38) in the form 


| gtg^gjaoL 

Jo 


where 



(39) 


(40) 


is regarded as the unknown function. Then equation (39) will he 
solved if any choice of g(o 1 ) can he found that will cause the 
integral I to vanish. 


The transformation 




g(a 1 ) = cn 2 h(e) 


a x 2 = n 2 sin 2 0 


a 2 = 


1 2 -1 


J 


converts equation (39) to 


I = C(l-q 2 ) 2 


r 

Jc 


Jt/2 


h(0) sin 0 cos 0 d0 
(sin 2 0 + q 2 cos 2 0) 2 


(41) 


(42) 


Certain s imil arities to integrals 14 and 15 of reference 11, table 48, 
suggest 


h(0) = tan 0 
h(0) = cot 0 


(43a) 

(43b) 
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as possible solutions of equation (38) or (39) . Both functions 
are found upon substitution and integration to satisfy the 
equation. Any linear combination of solutions (43) is also 
clearly a solution. Thus, the integral equation (38) does not 
have a unique solution. 

The later identification of equation (43a) as the physically 
correct solution is anticipated at this point. Then the u -velocity 
in the w = 0 regions of flow in may be evaluated as follows. 
Equation (36) can be written 


c^u _ + 4xg ^ 

^ rtp/^l-C^ 


(44) 


where I is the integral defined in equation (39) or equation (42) 
With the solution equation (43a) substituted in equation (42), the 
integral I has the value (reference 11, table 48, integral 14) 


I = C(l-q 2 ) 2 JL 


n < |a| < 1 


But according to equations (41) 




whence 


I 


rt 

4 


Cn 4 

o(o2- n 2) 3/2 


n < |a| < 1 


and, upon returning 'to equation (44) , 




4 

Cn x 

Ha 2 -n 2 ) S/S, IZ7 


(45) 


n< fa| < 1 


(46) 
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The u velocity in the region n < a £ 1 results upon inte- 
grating equation (46) twic6 between the limits O and 1. This 
integration may he effected with the aid of the elliptic -integral 
substitutions of appendix C, and the result is 


u(q,x) = 


Cn 2 x 

P(l-n 2 ) 


E(cp,k)-n 2 F(cp,k)j 



sin 2cp\ 


where 


cp = sin' 


-1 


/ l-O 2 

1-n 2 


k 



The value of C is established by the condition that 
u(n,x) = Ky = Knx/6. The negative range of cr (-1 <_ ct <_ -n) is 
taken care of by introduction of the antisymmetry factor a/|a| . 
These operations yield the final result 


n < |a| < 1 


(47) 

Equation (47) is the solution for the u velocity in the w = 0 
sectors of flow HI (fig. 2). This equation results from the choice 
of equation (43a) among the solutions of the integral equation (39) . 
The correctness of this choice is determined by comparison of 
equation (47) with the solution equation (15) of the related flow HI 
of figure 2. Both equations exhibit a half -order singularity in 
du/Sa at a = n+ together with du/da =0 at o = 1. 

(Equation (43b), on the other hand, leads to a different behavior.) 
Equations (47) and (15) are platted together in figure 6 for 
n = 0.707. Again, as for flews I and I* and. H and H * the agreement 
is such that the two sets of points lie on a single curve, the 
bottom curve in the figure. The degree of agreement as a function ' 
of the par ame ter n for the ray defined by cp = 60° is shown in fig- 
ure 7. The error is within 5 percent down to n = 0.4. 


u(ct,x) = 


I a l P 


|a| |E(cp,k)-n 2 P(cp,k)] 




sin 2Cp 


E 


(i’ k ) -“Hf ‘ 
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Flow IT 

The boundary conditions (fig. 1) show that this flow is quasi - 
conical like flow 111 . That is, the velocity components are of the 
foim u = xF(y/x) and w = zG-(y/x). The flow is, however, symmetric 
like flow X. A flow of this type can "be built up (ccmpare, treat- 
ment of flows I and 111) by superposing elementary source -sheet 
sectors that have a strength variation proportio nal to x. These 
sheets are arranged in the same manner as the correspo ndi ng constant - 
strength sheets employed for flow I. (See section A-A, fig. 5.) 

Let the source strength of an elementary sheet sector be 
specified by 

. w = (Af) x 

where Af is a constant for a particular sheet and the sector 
extends from a = - ap to o=Oj_. The incremental u velocity 

induced by such a sheet has been calculated by the method of 
reference 5 and the result is 

Au = WO^a,^) (48) 


where 






cosh"-*- 


1-OjO 

I a l- a l 


l+c^a 

Vl -° i 2 


cosh' 


-1 


l+c^a' 
a l +a , 


nx 




l-a^a 

-_L - - + cosh x 

°l“ a 


1+010 
0!+al 


I dx 


Mach 

cone 
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(The Integral -was not evaluated, inasmuch as only the z -derivatives 
of ¥ Trill he needed in the applications herein.) 

The individual sheet-strength parameter Af can he expressed 
as 


Af - ' d§7 to i 


so that 


df 

dCFj_ 

Au = - — dC^ 


The total value of u contributed hy the superposition of all the 
contemplated sheets is then 


rn 


u(x,a) 


1 

P 


JO 


df 


a^-w( x,^) a0 x 


(49) 


The development beyond this point runs similarly to that for 
flow IU. Inasmuch as u(x,cr) is specified over the shaded area 
(|al< n) of flow IV of figure 2, equation (49) is an integral 
equation for the unknown function df (C7^)/dCJ-^. The kernel 

¥(x,CJ,a;]_) is disagreeably complex, but its second derivative 

with respect to x. 


^■w(x,a,oi) 
ox 


-4d 4 a 1 


j txta^-c^) 2 


is relatively simple. The theorem, equation (D4), for differenti- 
ation under the integral sign is therefore employed to give 


\2 \ 

a u 


1 

F 


n 


i J£ a. 


(50) 
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The integration symbol r; has the special significance previousls 

Jo 


discussed under flow HI. 


The prescribed condition on u in the range |a| < n is 
(flow IV, fig. 2) 




U = KX 


or 


. CK 2 


> |CT| < n 


J 


Application of this condition to equation (50) yields the integral 
equation for df(o^)/d 0 ^ 




r n _ar 
i fx 5 

Jo (Ol 2 -^)' 


I cr | < n 


(51) 


This equation may be compared with the corresponding equation 
for flow 111 (equation (38) ) . Both equations are of the f orm of 
equation (39), but with different expressions for the unknown 
g(CT^). Thus the method of solution for equation (39) will apply 

for equation (51), with 


■ af («> 

O 

Once again the proper solution for g(a^) - Cn h(0) is 

found (by later comparison of flows IV and IV’) to be given by 
equation (43a) . This solution is now employed in equation (50) 

in the range n < |a| <1 for which c^u/ck 2 does not vanish. 

The integral has the value given in equation (45), and equation 
(50) becomes 
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gu 

,dx 2 


Cn 4 a 3 


y Px(a 2 -n 2 ) 


372 


aJi-o 2 


n< |a| <1 


(53) 


Now u is known to "be of the form (quasi -conical flow) 


so that 


u = Pyf (a) 


= - c^f ' (a) = g(a), say } 


d 2 u\ a ( / \ 


(54) 


Thus "by comparison of equation (53) with the last of equations (54) 


g* (a) 



The elliptic -function substitutions of appendix C axe convenient 
for performing the integration for g(o) , with the limits O 
and 1. Then f*(a) = -g(a)/j 2 may in turn he integrated to 
yield u. The result is 


u(CT,x) 


|F(cp,k)-E(cp,k) 
P(l-n 2 ) L 


where 


cp = 
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The constant C is established by the condition that 

u(n,x) = kx 


so that finally 


u(a,*) .uxg foH 

F 


2 ’ k ) “ E P k 


(55) 


Equation (55) is the solution for the u velocity in the 
w ss 0 sectors of flow XV of figure 2. This equation results 
frcm the choice of equation (43a) among the solutions of the 
integral equation (39) . The correctness of this choice is deter- 
mined by comparison of equation (55) with the solution equation (17 ) 
of the related flow 17* of figure 2. Both equations exhibit a 
half -order singularity in du/da at a = n+ together with 
du/dcr » 0 at 0=1. (Equation (43b), on the other hand, leads 
to a different behavior.) Equations (55) and (17) are plotted 
together in figure 6 for n = 0.707. Once again, aB for the pairs 
of flows I and I', II and II * , and III and HI*, the agreement is 
sufficiently close so that the two sets of points seem to define a 
single curve, which is the curve next to the bottom in the' figure. 
The degree of agreement as a function of the parameter n, for the 
point defined by cp = 60°, is shown in figure 7. The error is 
within 5 percent down to n = 0.27. 


DISCUSSION 

In the foregoing derivations principal attention has been 
devoted to the solution for the u velocity (or the pressure dis- 
turbance) in the regions in which it is unknown. The solution for 
the w velocity (or the upwash) in the regions in which it is 
unknown has also been obtained, although for brevity it has been 
emitted herein (except in the single example, appendix B). Com- 
parisons of the singularities in the w velocity between flows I 
and I*, HI and HI’, and 17 and 17*, respectively, were made in 
add i tion to the u velocity comparisons mentioned earlier. These 
additional comparisons again showed very close agreement between 
each flow and its related flow, in support of the correctness of 
the choice in each case among the solutions of the integral 
equations for flows I, IH, and 17. 
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A primary aim of this report has "been to provide the analytical 
basis for approximating in simple fashion the effect of subsonic 
trailing edges on damping in pitch and roll for thin sweptback wings 
in a supersonic stream. (See reference 3 for detailed application.) 
Accordingly, the main emphasis has been placed on the "desired flows " 

I to IV. - The more easily solved "related flows" I* to IV* were found, 
however, to approximate flows I to IV in the regions of interest 
much better than had been anticipated. (See figs. 6 and 7.) It is 
clear from the comparison that a moderate modification' of the boundary 
condition near the left-hand Mach line (fig. 2) - the distinguishing 
feature between the unprimed and primed flows - will scarcely affect 
the flow in the region near the right-hand Mach line. Thus such a 
modification may freely be made to simplify a given problem, and the 
solution of the modified flow will apply to the desired flow with 
engineering accuracy . In particular, the related flows herein may 
be used in place of the desired flows if trigonometric functions are 
preferred to elliptic integrals. 

A general method for the calculation of flows (fig. 9) of 
which flows I* and IV* are special cases is Included in a paper by 
Goodman (reference 12) after the bulk of the present work was 
completed. The method and its applications have been simplified 
and extended by Mirels (reference 13) . The general solution for 
the unknown u velocity in the flows specified in figure 9 is given 
In elegant form in reference 13. This result, originally obtained 
by means of a doublet distribution, may be derived by means of a 
development of the present source -distribution method. The unknown 
w velocity may be obtained as well. The details are given in 
appendix E as a matter of interest. 


COECLUDHTG REMARKS 

The flow over a sveptback wing may be obtained by superposing 
on a basic delta-wing flow additional, flows to cancel the lift 
outside the boundaries defining the sweptback wing. The cancel- 
lation flow for the trailing edge modifies the lift in a region 
ahead of the trailing edge if the component of the stream velocity 
normal to the edge is subsonic. For angle of attack the principal 
part of this cancellation flow is conical; for rolling and for 
pitching the principal part is quasi -conical. The derivation of 
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these and several related conical and quasi-conical flows has 
been carried out in this report. For each case the problem was 
formulated as an integral equation. Some of the results are 
applied to damping in roll and in pitch in another report . 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 3, 1950. 
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APEEKDIZ A 
SYMBOLS 


The following symbols are used in this report: 

C houndary for line integral, as specified in text; also 

undetermined constant 

0^ pressure coefficient (proportional to u) 

E(cp,k) incomplete elliptic integral of second kind with amplitude 
cp and modulus k 

F(cp,k) incomplete elliptic integral of first kind with amplitude 
<p and modulus k 

f(0) function defined in equation (22) 

f (a) function defined in equation (54) 

f(a^) function defined in equation (33) 

g (a) function defined in equation (54) 

g(a j) function defined in equation (40) for flow XU, in 

equation (52) for flow XV 

h(0) function defined in equation (41) 

I integral defined in equation (39) 

J integral defined in equation (Bl) 

k = \jl-n 

M Mach number, ratio of stream velocity to velocity of 

sound in free stream 

n value of O along side edges of shaded triangular regions 

in figure 2 

q = a/^O 2 - n 2 

S area -for surface integral, as specified in text 
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uC^o-l) 

U,v,v 

u 0 

Au 

V 

W(x,a,oi) 

X,J,Z 

a 

P 

5 

e 

k 

a = 
t = 

cp = 

i,n 

ih,iy 

iSu*, 

HI',IY' 


function defined following equation (34) 

disturbance velocity components along x-, y-, and 
z-axes, respectively 

constant value of u over prescribed area 

increment in u 

free -stream velocity 

function defined following equation (48) 

Cartesian coordinates: x-axis parallel to free -stream 

direction; y-axis horizontal and toward right, looking 
upstream; z-axis vertically upward 

angle of attack 

Vm 2 -i 

defined in equation. (22) 

defined in equation (22) for flows I and II, in 
equation (41) for flows HI and IT 

constant of proportionality (fig. 2) 

oblique coordinates measured parallel to downstream Mach 
lines, defined by equations (2) 

P y/x 

j— (See fig. 3 for geometric significance.) 



desired flows designated in figure 2 


related flows designated in figure 2 
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n> 


special integration sign defined following equation (36) 


A primed function signifies the first derivative of the 
function with respect to the independent variable. 
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APPENDIX B 


EVALUATION OF UPWASH TO BEGION B OF FLOW I’ 

Note that Just "before equation (4) there is the definition 


“ ii 


u 


V 7 


VRT 


(Bl) 


Now J is a known function, which was evaluated in the course of 
obtaining equation (10) . Thus equation (Bl) is an Abel integral 
equation for the unknown function dw/dx^. Its solution (refer- 
ence 6) is j 


n 


i 


^W _ m d 

^ 311 


V T 




(B2) 


With the value of J for flow I* (equation (6)) this 1 b, provided 
that the integral on the right-hand side is continuous (reference 7) 


dw 


MUq 

* 



V5=l 


^0 


Upon dropping the subscripts and converting to Cartesian coordinates 
by means of equations (2), 


dw = _ Pu 0. ) 2(l-m) 

it ^(x+Py)(nz-Py) 


(B3) 
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Tlie indefinite integral is 


v = 


^0 ) 2 (1-t-n) cosll -l 2mc-py(l-n) 
]j n f3y (1+n) 


fed+n) 
V n 


cosh 


-1 


2n-g(l-n) 
|a| (1+n) 


(B4) 


Equation (B4) gives v = 0 along the £-axis (a = -l) and along 
the line q = t£ (a = n), as it should for flow I'. The constant 
of integration is therefore zero and equation (B4) is the required 
solution far the upwash velocity v. 


1317 



1317 


NACA T3J 2147 


37 


APEEKDK C 


AIDS TO IUTEGEATION 

The more elaborate Integrations in this report, leading in 
most cases to elliptic integrals, are most easily evaluated with 
the aid of the following elliptic -function substitutions, together 
with the tables of integrals in references 14 and 15: 



= k sn u 


i\jo 2 -n 2 


k cn u 


a = dn u 
2 

da = -k sn u cn u du 


where 


u a F(<P.,k) [not to be confused with the 
u velocity] 

k = ^ 1-n 2 
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APPEUDIX D 


OF HJTEGBAIS OF SINGULAR FUNCTIONS 


By Franklin E. Moore 

The problem arising in "the differentiation of equation. (35) 
may he stated in the following general terms : It is required to 

differentiate, -with respect to a parameter of the integrand, the 
definite integral of a function that, though integrable, has a 
singularity inside the range of integration suoh that the definite 
integral of the derivative is not convergent. It is shown herein 
that, subject to oertain restrictions, this differentiation can be 
carried out, in a simple manner, avoiding consideration of the 
singularity. 

The integral to be differentiated can be written 

"b 

I(r) = f(x,!)d£ 

Ja 


It will be supposed that f (x, |) has a singularity at some point 
a < | < b and is a function such that 


(1) An indefinite integral F(x,t) = 


>1 


f(x, !)d£ may be 


found. 


(2) The function l(x) is a convergent improper integral that 
can be written l(x) = F(x,b)-F(x,a) . (See reference 16, para- 
graph 169.) Under these restrictions, f(x) may have an "integrable 
singularity" within the interval of integration. For example, 
f(x,|) = logg |x- 1 1, a < x < b, meets the foregoing requirements 

(1) and (2). 


It is required to carry out the following differentiation: • 
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By definition. 


fix, t) 




F(x,t) 


It is supposed that 


(3) At the points £ « a,h, F(x,£)J =» 

In order that this relation he valid, the following restrictions 
on F(x,l) in the neighborhoods of £ = a,b are required (ref- 
erence 16, paragraph 213) : 


(a) 


and exist. 


(b) «iL * s continuous, 


It then follows that 


^ f(x,t,) a 




Thus, provided that 


(4) ^rf(x,|) has an indefinite integral. 





^ f(x,|)i + g(x) 


where g(x) is an arbitrary function. 
Equation (Dl) can then be written 


I*(x) = 


rt 


b 





f (x,|)d|+g(x) 

- 


— 

l-b 



^ f(x,l)d|+g(x) 


— * £=a 
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or 


i*fc0 = 


pa 


£f(x,g)d| - 


SE f(3C ' g)ag 


When the convention introduced in equation (36) is followed, 
equation (D2) can "be written 


(D2) 


' I’Cx) 





(D3) 


The result (D3) may he stated as follows : Subject to the 

restrictions (1) to (4) on the behavior of the integrand near 
the limits, the differentiation with respect to a parameter to 
the integrand of a convergent improper integral may be accomplished 
by formally integrating the derivative of the integrand, as though 
its singularity were not present. 

It is clear that 

r 

l"(x) = l if(x,t)it (D4) 

jja * 


provided that requirement (3) is replaced by the restriction that, 
near £ = a,b 



If 


a 

3T 


Ji 

Sbc 2 


F(x,S) 


Equation (D4) 


is the relation used in this report. 


The foregoing proof may easily be extended to the case where 
a and b are functions of x, subject only to the additional 
requirement that 
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(5) a'(x) and VCx) are continuous, yielding the result 


n 

I*(x) = h , (2:)f(x,'fa)-a'(3:)f(2:,a) + j Jr f(x,£)d£ 

J a 


(D5) 
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APPENDIX E 


GOOIMAN-MIKELS GENERAL IZE D CANCELLATION FLOW 


The generalized cancellation flow introduced in reference 12 
and extended in reference 13 is specified "by the boundary condi- 
tions of figures 9(a) or 9(1)* depending on the position of the 
line Ti = f(£). Consider first figure 9(a). The flow can he. 
represented, as before, by sources distributed along the surface 
in proportion to w. In the Mach coordinates of reference 5, the 
surface potential is 


qp( 1^) 



region B (El) 


where (|,tj) is located at P in region B. In what follows <T> 
will play essentially the role that u played in the solution 
of flows I* to IT' . 

Consider next the case of figure 9(b) . By virtue of the left- 
hand u = 0 region (reference 5), the region of integration for 
cp( £,t]) will- not include the area left unshaded in the upper comer. 
It is found that with the present notation the limits of integration 
axe unchanged from the case of figure 9(a) . Thus equation (El) 
applies equally well for the case of figure 9(a) or 9(b). 


Solution for u in region A . - Define 

n t 




U 




wdl-L 

VHI 


(E2) 


so that 


'P(L’l) = “ 


i f(!) 


V^i 


region B 


(E3) 
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Equation (E3) is an Abel integral equation for J. Its solution 
is (references 6 and 7) 


HUrti) - - - 


1 5 
jt 3 t £ 




f(£) 


cp(£,z)dz 

^l" 2 


(E4) 


Wow let the general point (£,q) be located in region A 
at Q. The only change in the expression for <P is in the upper 
limit for q^. This upper limit is now at the intersection of 

the line PQ with the line q = g(£). Accordingly, 


^g(£) 

cpCI^-n) = - 

Jf(t) 


j(|,T) 1 )dT ] 1 


(E5) 


The value of J already obtained in equation (E4) applies here 
as well as earlier. Upon making the substitution. 


cp(|,Ti) 



^g(t) 


■^1 


1 

rt 

dq-L 

V^l 

a 


cpd.zidz 

\j^l~ z 


f(|) 

vj 

f(l) 



region A 


(E6) 


P6«) 


ptli 



^1 


f (£) 


cp( £,z)dz 
Oil-*) 5/8 


region A (E7) 


where the sign | designates the finite part of the integral 

according to Hadamard's definition (reference 17). 
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The factor l/ which does not contain z, may he 

included in a common integrand, and the finite -part sign may he 
removed to the outside of the double integral: 


CP(M) - - i 


'g(t) 


dTh 


f(8) 


'^1 


f(|) 


c P(^,z)dz 




^72 


(EQ) 


The order of integration may now he reversed with the aid of the 
scheme for the limits sketched in figure 10: 


<pU,u) = - 


1_ 

2rt 


Pg(l) 


f(8) 


Cp( £,z)d z 


'g(l) 


dTli 


/\/ -n-ni (tii-z ) 3 / 2 


(E9) 


Evaluation of the finite part of the inner integral yieldB the 
final result 




'g(l) 


cp(|,z)dz 


f(8) 


(t]-z) \Jg(£)-z 


region A (E10) 


Equation (E10) gives the surface potential in region A that results 
from a prescribed distribution of surface potential in region B. 

Equation (E10) is equivalent to equation (15) of reference 13, 
which applies to the mirror image of figure 9. (The same equation 
was obtained for a restricted situation in reference 18 (equation (61d) 
therein) . ) According to reference 13, differentiation yields 
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u(e,i) 


= Vti-s(S) 


Pad) 




f(!) 


u(£,z)dz 

(n-z) \js(V~ z 


region A 


(Ell) 


which is equation (17b) therein. This equation constitutes a 
general solution for the disturbance in region A regardless of 
whether the Kutta condition is imposed along the line q = g(|). 

It is shown in reference 13 that imposition of the Kutta condition 
along the line q = g(|) modifies the sidewash distribution v 
in such a way as to cause the second integral to vanish. (If the 
boundary = g(£) of fig. 9 slopes everywhere toward the left, 
instead of toward the right as shown, v is determined uniquely 
by the known u distribution, and the Kutta condition may not 
be imposed.) 

The Kutta condition version of equation (Ell) (that is, with 
the second integral set equal to zero) may be obtained more directly 
for case (b) of figure 9. In this oase (including the special 
case f(|) = 0), u may be written in place of cp and dw/ckj 

in place of w in equations (El) through (E10); the result is 
given in the new equation (E10) . This approach at first parallels 
and then extends that of equation (3) and the following equations 
in the text. 

Solution for w in region B «- It Is worthwhile to complete 
the solution of the Goodman -Mirels cancellation flow (fig. 9) by 
evaluating the upwash velocity w In region B . This knowledge 
of w wl.11 be useful in determining the downwash field in the 
general vioinity of the wing - in particular, within a chord length 
behind the trailing edge - by means of the cancellation technique. 

(A line vortex method (references 19 and 20) yields the downwash 
farther back more simply, and simple means are known (cited in 
references 19 and 20) for determing the downwash immediately behind 
the trailing edge.) 


l-g'(£) 

2 JtP 




'g(|) 


f(t) 


Pu + V 
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dz 
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The solution is started "by inverting (solving) the Abel 
integral equation (E2) with the aid of references 6 and 7: 


v = M 




\ 




j(^T) 1 )dl 

WT 


(E12a) 




j(£,iii)at 


(E12) 


The solution for J has already "been determined in equation (E4) . 
Changing z to tj therein gives 




i a 

* ^7 


p^i 


f(6) 


qp(| ,-n)d.Ti 


(E13a) 



P^l 


f(l) 


gpC E, Tl)dTl 

(%-n) 3 ^ 


(E13) 


Substitution of equations (E13a) and (E13) into equations (E12a) 
and (E12 ) , respectively, yields 
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region B 


(E14a) 


4rt 


'Si 


'&(%) 


'•l 


— AL„ ... 

(5i-e) 3/2 






region B 


(E14) 


By comparison frith equation (10) of reference 13, equation 
(E14) may he interpreted as the upwash velocity w due to a 
given distribution of doublets (or of vortioity) over a prescribed 
area. The area of integration is shown cross-hatched in figure 11 
for the same two configurations (a) and (b) considered in figure 9. 
Note that the regions A' and B', which lie in the zone that influences 
the point are not included. This result implies that for 

the point of the doublets in A' exactly cancels 

the integral of the doublets in B'. The situation here with respect 
to doublets presents an interesting parallel to E-ward' s original 
discovery with respect to sources in connection with his membrane 
concept. (See reference 21.) 

Equation (E14a) is an alternate form of equation (E14) t hat 
may yield simpler integrations in some cases. Note that the [ 

(finite part) operation is avoided. Similar considerations apply 
to the alternate forms (E12a) and (E12), (El3a) and (E13) . 


A particularly simple result for ctr/dx can be obtained 
when the Eutta condition is imposed. Thus, apply to equation (E14) 
a procedure analogous to that by whioh dqydx is obtained from 
cp in reference 5. The result is 




1317 


NACA TN 2147 


49 


The line integral may be put in more convenient form by means 
of an integration by parts. The final result may be -written 


_ _ _M_ 

4rt 
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Jg-Ct)-, ) 


d| 

a 1 -i ) 3/2 


r r h 


if(£) 


u. ( ^ . TI ) d- 
(Tli-ti) 
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3/2 


1 - 




dm 




ffta) 


Pu+v 

\l8(i 2 )-r] 


cLri 


(E16) 


The line integral is the same as the one in equation (EU) . Here 
again, then, the line integral vanishes if the Kutta condition is 
imposed along the line £ = G-(t])~ 
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Figure 1. - Sweptback wing showing basic delta-wing load 

distribution along section A-A and approximate cancellation 
of this load behind wing by superposition of special flows. 
Special flows are plotted with reversed sign. Cross- 
hatched areas represent induced changes in loading ahead 
of trailing edge. 
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III III 



Figure 2. - Flows used for cancellation of lift behind trailing 
edge of sweptback wing, I, II, III, and IV, and related flows 
I', II', III', and IV'. 
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Figure 4. 


Data for flow II* • 
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Figure 5. - Superposition of elementary source-sheet 
sectors to build up flow I. (Schematic) 





Figure 6. - Comparison of distribution of u velooity (proportional to lift) in right-hand 
vr=0 regions (n < o < 1) of flows of figure 2. n=0.707. 
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Velocity ratio 






Figure 8. - Superposition of elementary source-sheet 
sectors to build up flow III. (Schematic. ) 













